Gaussian expansion approach to Coulomb breakup 
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An accurate treatment of Coulomb breakup reactions is presented by using both the Gaussian expansion 
method and the method of continuum discretized coupled channels. As L 2 -type basis functions for describing 
bound- and continuum- states of a projectile, we take complex-range Gaussian functions, which form in good 
approximation a complete set in a large configuration space being important for Coulomb-breakup processes. 
Accuracy of the method is tested quantitatively for 8 B+ 58 Ni scattering at 25.8 MeV. 



PACS numbers: 24.10.Eq, 25.60.Gc, 25.70.De, 26.65.+t 

Determination of the neutrino oscillation parameters is one 
of the central issues in the neutrino physics. The astrophysical 
factor 5i 7, which is essentially equal to the cross section of the 
p-capture reaction 7 Be(p, 7) 8 B, plays an important role in the 
parameter-search procedure, since the prediction value for the 
flux of 8 B solar neutrino is proportional to S17. The required 
accuracy from astrophysics is about 5% in errors 1 1]. Due to 
difficulties of direct measurements of 7 Be(p, 7) 8 B at very low 
energies, alternative indirect measurements were proposed; 

8 B Coulomb breakup QUHBHUisa typical exam P le 
of them. 

The method of continuum discretized coupled channels |8, 
9] (CDCC), which was shown to describe various proiectile- 
breakup processes 

G3 EH 112 IS HIS IE El, 

has 

also been applied to 8 B Coulomb breakup with high degrees 
of success 1 7, 19, In the all analyses above, CDCC 

describes the reaction system with a three-body model, i.e., 
p+ 7 Be+A as shown in Fig. 1. 
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FIG. 1: Illustration of a three-body (A+b+c) system. The symbol 
B=b+c stands for the projectile and A is the target. 

From 8 B Coulomb breakup measured at intermediate en- 
ergies, S17 was determined using the first-order perturbation 
theory QllBlllllI It was shown in Refs. HHIIl that S17 
can alternatively be determined by CDCC analysis, combined 
with the asymptotic normalization coefficient method 112 311 . of 
8 B breakup even at low energies, such as the Notre Dame 
(ND) experiment ll24ll . It was found, however, that the ND 
data contain an irrelevant component to S17, corresponding to 
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the l/2~ 7 Be excited state in the final channel, which of or- 
der 10% presents in the 8 B ground state 12511 : effects of this 
component on Si 7 can be expected around 10%. Therefore, 
the evaluated value of S17 in Refs. l2"lH22ll should be assumed 
as the upper limit of S17. Even though the irrelevant compo- 
nent can be experimentally removed by detecting the emitted 
7 Be in the ground state, the 1/2 - 7 Be excited state may play 
a dynamical role during the breakup process of 8 B. Thus, de- 
scription of 8 B breakup beyond the three-body model, i.e., 
p+ 3 He+ 4 He+A four-body model, is highly expected in order 
to determine Si 7 with very high accuracy. 

Very recently, a new treatment of breakup continuum in 
CDCC was proposed Eo], which uses pseudo-state (PS) 
wave functions 0, 0, l28ll obtained by diagonalizing inter- 
nal Hamiltonian of the projectile with Gaussian basis func- 
tions 1 29]. For nuclear breakup processes, CDCC with the PS 
method (PS-CDCC) was found to perfectly reproduce the "ex- 
act" breakup S-matrix elements S(fc) calculated by the stan- 
dard CDCC, i.e., with the Average (Av) method llBlSIIl 
I3II1 for discretization of breakup continuum. The most impor- 
tant feature of PS-CDCC is that it generates S(k) as a smooth 
function of k in its entire range, without assuming any form 
a priori for /c-distributions, in contrast to CDCC with the Av 
method (Av-CDCC). 

Additionally, as discussed in Refs. jH HU, PS-CDCC 
with Gaussian basis functions makes it possible the four- 
body CDCC analysis of projectile breakup, a preliminary re- 
sult of which for 6 He elastic scattering, including effects of 
breakup channels of 6 He, has already been obtained 1321 . Be- 
fore making four-body PS-CDCC analysis of 8 B Coulomb 
breakup, however, one must see the applicability of PS-CDCC 
to breakup processes due to long-ranged Coulomb coupling 
potentials. This is quite nontrivial because the r-space of the 
projectile needed to describe Coulomb breakup is much larger 
than that in the case of nuclear breakup; for example, the ef- 
fective space is r < 20 fm for nuclear breakup of dor 6 Li 1 26], 
but r < 100 fm for Coulomb breakup of 8 B IHEl . 

The aim of this brief report is, as the first step towards the 
four-body CDCC analysis of 8 B Coulomb breakup, to show 
that three-body PS-CDCC based on the Gaussian basis func- 
tions can well describe the corresponding result of three-body 
Av-CDCC for 8 B breakup. As for the test case we take 8 B 
breakup from 58 Ni at 25.8 MeV. 
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Below we briefly recapitulate the formulation of three-body 
CDCC. We consider the three-body system of Fig. 1. The 
model Hamiltonian of the system is 

H = K r + V hc (r) + K R + U hA (r hA ) + U cA (r cA ) ■ 

Coordinates are defined in Fig. 1. Operators K r and Kr are 
kinetic energies associated with r and R, respectively, and 
Vbc(r) is the interaction between b and c. The interaction Ub A 
(U cA ) between b (c) and A is taken to be the optical potential 
for b+A (c+A) scattering; Coulomb breakup is induced by the 
Coulomb components of the potentials. For simplicity, in this 
study we neglect the intrinsic spins of individual constituents 
of the system. 

In CDCC, the states of the projectile are classified by the 
linear and the angular momenta, k and £, of relative motion 
between b and c, which are truncated by k < fc max and 
£ < ^max- The truncation is the most basic assumption in 
CDCC, and it is confirmed that calculated 5-matrix elements 
converge for sufficiently large fc max and £ max Mil (ilSOIl. The 
converged CDCC solution is the unperturbed solution of the 
distorted Faddeev equations, and corrections to the solution 
are negligible within the region of space in which the reaction 
takes place [33]. The fc-continuum of the b+c system are then 
discretized into a finite number of states, each of which corre- 
sponds to a "discretized-continuum state" with a certain posi- 
tive eigenenergy labeled by i. The resulting orthonormalized 
wave functions, {&u(r)i e Ye m (£l r ); i — l—N}, are assumed 
to form an approximate complete set in the r- and fc-space 
being important for the breakup reaction. 

The three-body wave function ^ of the system is expanded 
by {&i£(r)i Yg m (Q, r )} and then inserted into the approxi- 
mate three-body Schrodinger equation (H — E)^ = 0. One 
can then obtain a set of coupled differential equations, called 
CDCC equations, that provide the discrete 5-matrix element, 
5 7 , for the transition from the initial channel to a discretized- 
continuum one 7 = (i,£). In order to calculate a breakup 
cross section of the projectile with a certain range of breakup 
energies, or a coincidence-cross-section, a smoothing proce- 
dure, which constructs continuous Se(k) from the discrete 5 7 , 
is necessary. Actually, this is the case in the analysis of 8 B 
Coulomb-breakup experiments 0, H . 

Discretization of the breakup continuum in the PS method 
is done by diagonalizing the internal Hamiltonian iJbc = 
K r + Vbc(r) in a space spanned by a finite number of i 2 -type 
basis-functions, for which we here take the following pairs of 
functions [29], 

4>%{r) = /exp [-(r/a,) 2 ] cos [ b (r/a,) 2 ] , 

4>%(r) = /exp [-(r/a,) 2 ] sin [ b (r/a 3 ) 2 ] , (j = 1-n), 

where {a 3 } are assumed to increase in a geometric progres- 
sion and b = tt/2. We refer to the basis as the complex- 
range Gaussian basis, since the basis functions can be ex- 
pressed by Gaussian functions with a complex-range parame- 
ter, r l exp[— (1 + ib)(r/dj) 2 ], and its complex conjugate. The 
complex-range Gaussian basis functions are oscillating with r, 
so they can simulate the oscillating pattern of the continuous 



breakup-state wave-functions as shown in Ref. l29Tl . which is 
very important for the description of Coulomb breakup by PS- 
CDCC. An accurate transformation from fiL to Si (k): 

StW^J^fuik)^, (1) 

i 

with fu(k) = (fc, r)\&ii (r)), is possible when the basis 
functions form an approximate complete set in the finite con- 
figuration space being important for the breakup process 1 26]. 

In the Av method, on the other hand, the fc-continuum [0, 
fcmax], for each £, is divided into a finite number of bins, each 
with a width An = ki ~ and the continuum breakup- 
states in the i th bin are averaged with a weight function 
wu(k) |]8j, IPj ; we in this study take wu(k) = 1 for simplicity. 
The resulting orthonormal state is then described as 

$ u (r) = -£= ^ Mk,r)dk (forAv). (2) 

V &U Jki-! 

Inserting Eq. (|2} into Eq. Q leads to S^ v (k) — S ina /y/Au 
for fci_i < fc < ki, as shown in Ref. II 911. 

It should be noted that Sf LV (k) is smooth only within each 
bin, while Si(k) with the PS method, i.e., S^^fe), is a smooth 
function of k in its entire range. Details of the discretization 
and smoothing procedures both for the Av and PS methods 
can be found in Ref. [26]. 

In order to see the applicability of PS-CDCC for Coulomb 
breakup, we calculated the breakup cross section of 8 B+ 58 Ni 
scattering at 25.8 MeV. In this process both nuclear and 
Coulomb interactions, with interference, play important roles. 
Moreover, first order perturbation theory cannot reproduce 
the experimental data [24] at all, since higher-order processes 
cannot be neglected at such a low incident energy 1 34] . Thus, 
it is a good test for PS-CDCC to compare the results with 
those calculated by Av-CDCC that reproduce the experimen- 
tal data very well (13 EH Q . 

For both the Av and PS methods, we took the 8 B single- 
particle wave-function of Esbensen and Bertsch |35] and the 
same potential as for the ground state was used to generate 
wave functions in scattering states; the depth of the potential 
is fitted to reproduce the separation energy of the proton, i.e., 
137 keV. We included only s- and p-states of 8 B for saving 
computation time. It was shown in Refs. flDl EH l22ll that 
d- and f-states are necessary to reproduce the experimental 
data 12411 . Our aim in this brief report, however, is to see 
the convergence between two theoretical calculations. As for 
the distorting potential forp- 58 Ni ( 7 Be- 58 Ni), we adopted the 
parameter set of Yamashita 12111 (Moroz et al. |36]). The 
maximum total-angular-momentum J max is 1000 and i? max 
is taken to be 500 fm. CDCC equations were solved with the 
predictor-corrector Numerov method with stabilization l37ll . 

In the Av method, the modelspace with k max = 0.66 fm -1 
and An = 0.66/16 (0.66/32) fm" 1 for I = 1 (0) gives con- 
vergence of the resulting total breakup cross section. The 
maximum internal coordinate r max was taken to be 100 fm. 
In order to obtain the correct asymptotic form of the Coulomb 
coupling potentials, we made the following approximation for 
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FIG. 2: Discretized momenta for B; the left (right) side corresponds 
to the s-state (p-state). The horizontal dotted line represents the cut- 
off momentum femax taken to be 0.66 fm -1 . 
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good convergence. The number of channels included in the 
CDCC calculation, i.e., k < fc max = 0.66 fm -1 , was 18 for 
both the s- and p-states. The resulting wave functions with 
positive eigenenergies turned out to oscillate up to about 100 
fm. In order to obtain orthonormality of the $^(r), we put 
r ro ax — 130 fm in the calculation of coupling potentials. We 
show the level sequences of the resulting discrete eigenstates 
in Fig. 2. One sees each interval of the levels is almost even, 
which is just the same as in Ref. l2cill . 

We show in Fig. 3 the calculated angular distribution of 8 B 
total breakup cross section for 58 Ni( 8 B, 8 B*) at 25.8 MeV. The 
solid and dashed lines correspond to the Av and PS meth- 
ods, respectively. One sees both results coincide with very 
high accuracy, which implies that a forthcoming four-body 
PS-CDCC analysis of 8 B Coulomb breakup has enough ac- 
curacy to determine a reliable value of Sn, i.e., within 5% of 
errors. 

In order to see the validity of PS-CDCC for the 8 B 
Coulomb breakup more strictly, below we compare Sf s (k) 
with Sf^ v (k). In the calculation of the latter, i.e., the "exact" 
Se(k), we put An for the p-state to be 0.66/32 fm -1 . This 
makes the fc-dependence of S^ v (k) clearer, since S^ v (k) is 
constant in fc-region corresponding to each bin with the width 
of Am, as shown in Ref. |26]. It should be noted that this 
refinement of the modelspace for the Av method makes no 
changes in the physical quantities such as the elastic and total 



where j = p or 7 Be. In the practical calculation we took ac- 
count for Coulomb radii for p- 58 Ni and 7 Be- 58 Ni and Eq. (0 
has a slightly complicated form I21I1 . It should be noted that 
Eq. (|3j is exact when r max — > oo; the numerical results shown 
below were found to converge at r max = 100 fm. The result 
thus obtained, i.e., with Av-CDCC, is assumed to be the "ex- 
act" solution. 

In the PS method, we used the complex-range Gaussian ba- 
sis with (ai = 1.0, a n — 35.0, 2n = 60, b = ir/2) that gives 
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FIG. 3: Angular distribution of the total breakup cross section for 
58 Ni( 8 B, 8 B*) at 25.8 MeV. The solid and dashed lines represent the 
results with the Av and PS methods, respectively. 
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FIG. 4: The squared moduli of breakup S-matrix elements, as a func- 
tion of k, at J = 150 for s B+ 5S Ni scattering at 25.8 MeV. The 
panel (a), (b), (c) and (d) correspond to (£, L, Lo) = (1, 150, 150), 
(1,149,151), (1,151,151) and (0,150,151), respectively. In each 
panel, the solid line represents the result of PS-CDCC, while the step 
line is the result of Av-CDCC assumed as the "exact" S-matrix ele- 
ments. 
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breakup cross sections. 

Figure 4 shows the result of the comparison of |SV(fc)| 2 
at J = 150, which corresponds to the scattering angle of 
10° assuming the classical path. It was found that CDCC 
calculation with only Coulomb coupling potentials gives a 
peak at 10° in the total breakup cross section. Thus, it can 
be assumed that Fig. 4 corresponds to the most-Coulomb- 
like breakup process; in any case, features of the result were 
found to be almost independent of J. In each panel of Fig. 4 
the result with (I, L, L ), where L (L ) is the orbital angular 
momentum between 8 B and 58 Ni in the final (initial) chan- 
nel, is shown. The panel (a), (b), (c) and (d) correspond 
to [i,L,L ) = (1,150,150), (1,149,151), (1,151,151) and 
(0,150,151), respectively; all other components are negligi- 
bly small and not shown in the figure. One sees that the re- 
sult of PS-CDCC (solid line) very well reproduces the "ex- 
act" solution (step line) for all k being significant for the 8 B 
"Coulomb" breakup. 

In summary, PS-CDCC proposed in Ref. i26ll is shown to 
describe Coulomb breakup processes very well. Due to the 
long-ranged Coulomb coupling-potentials, the modelspace re- 
quired for CDCC is very large. Particularly, one must prepare 
the internal wave functions of the projectile, both in bound 
and continuum states, for a wide range of the internal coordi- 
nate, say, 0-100 fm, which is in general difficult for PS meth- 



ods. We found that this can easily be achieved by using the 
complex-range Gaussian basis, in the case of two-body pro- 
jectile. According to Ref. 1 29], the basis is also applicable to 
the reaction processes with three- and four-body projectiles, 
since energies and wave functions of the pseudo-states of the 
projectiles are given easily. Moreover, all coupled-channel po- 
tentials in four- and five-body PS-CDCC can analytically be 
given by expanding each nuclear-optical-potential concerned 
in terms of Gaussian functions. Then, PS-CDCC based on 
the complex-range Gaussian basis functions is expected to be- 
come an effective method of practical use even for four-body 
Coulomb breakup of 8 B, in which 8 B is assumed to be a bound 
state of the p+ 3 He+ 4 He system. Such a four-body analysis 
of 8 B dissociation is desirable in order to determine with 
high accuracy and then approach to the solar neutrino problem 
from the nuclear physics side. 
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